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This paper is dedicated to the memory of Ralph S. Phillips (1913- 
1998), one of the seminal figures in the theory of semi-groups. 
Early in our careers, he was a profound and positive influence 
on both of us. 

Abstract 

We analyze the diffusion process associated to equations of Wright-Fisher 
type in one spatial dimension. These are associated to the degenerate heat 
equation 

d t u = a(x)d x u + b(x)d x u (1) 

on the interval [0, 1], where a(x) > on the interior and vanishes simply at 
the endpoints, and b(x)d x is a vector field which is inward-pointing at both 
ends. We consider various aspects of this problem, motivated by their appli- 
cations in biology, including a comparison of the natural boundary conditions 
from the probabilistic and analytic points of view, a sharp regularity theory 
for the "zero flux" boundary conditions, as well as a derivation of the precise 
asymptotics of solutions of this equation, both as t — > 0, oo and as x — ► 0, 1. 
This is a precursor to our more complicated analysis of these same questions 
for Wright-Fisher type problems in higher dimensions. 

1 Introduction 

Consider the differential operator 

d 2 d 

L = a{x)— i + b{x)— (2) 
ax z dx 
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on the interval [A-B], where the coefficient functions a(x),b(x) G C°°([A, B]). 
Our main assumptions on the coefficients are that 



a{x) = (x — A)(B — x)a(x), 

where a(x) £ C°°([A, B]), and a(x) > OVx £[A,B] 



(3) 



and 



b(A) > 0, 6(5) < 0. 



(4) 



In other words, we assume that a(x) vanishes simply at the endpoints of the in- 
terval and nowhere else, and that the first order term b{x)d/dx is an inward pointing 
vector field. The diffusion associated to L is of importance in population biology. 
The basic example is the so-called Wright-Fisher operator 



on the interval [0, 1], which is the diffusion limit of a Markov chain modeling the 
frequency of a gene with 2 alleles, without mutation or selection. If the muta- 
tion rates between the alleles are /X12 and /421, and the mutant allele has selective 
advantage s, then 



see 0J. Accordingly, we shall call any operator of the form Q with coefficients 
satisfying (0) and (0]) a generalized Wright-Fisher operator. By an affine transfor- 
mation, we can also reduce to the case where L is defined on the interval [0,1]. 

In a seminal 1952 paper, Q, Feller considered the additional conditions one 
should impose on the domain of an operator like L to obtain the generator of a 
positivity preserving C° -semi-group. Feller's starting point was the Hille-Yosida 
theorem and his analysis centered on the construction of the resolvent kernel, (L — 
A) -1 . Our approach, by contrast, focuses directly on the Schwartz kernel for e tL . 
In addition, we restrict to a particular choice of boundary condition, which arises 
in applications to population genetics. 

The plan of this paper is as follows. In the next section we recall a change 
of variables introduced in [5j which reduces a general Wright-Fisher operator L 
to one with principal part x(l — x)d%. We then describe the natural boundary 
conditions for this operator and its adjoint. The operator L is modeled near each 
endpoints, x = or x = 1 by an operator of the form: 



L WF = x(l-x)— 5 



(5) 



b(x) = nuil - x) - \l 2 \X + sx(l - x) 



(6) 



Li = xdl + bd. 



(7) 
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where b is a nonnegative constant. The model operators act on functions on [0, do), 
with boundary conditions at x = induced from those for L. The next several sec- 
tions are devoted to a careful analysis of the solution operator for <9 t — L^. After a 
discussion in §|5]of maximum principles in this setting, we derive in §|6]an explicit 
formula for the Schwartz kernel k% of this solution operator, and then in §|7]com- 
bine these ingredients to prove sharp mapping properties for the kernels k\. The 
long §[8] contains a variety of technical estimates needed to analyze perturbations 
of these model operators. After these extensive preliminaries, it is straightforward 
to assemble this information and express the fundamental solution for a general 
Wright-Fisher operator L on [0, 1] as a convergent Neumann series. This is done 
in § |9j this leads directly to the description of precise asymptotics of solutions of 
(dt — L)u = as t — > and as x — > 0, 1. An asymptotic expansion of the heat 
kernel for L was obtained using other methods by Keller and Tier in [9]. The next 
section §[l0]uses this to provide a characterization of the infinitesimal generator of 
the induced semi-group on C° and its spectrum, which gives the long-time asymp- 
totics for solutions of the diffusion problem. In § [TT] we then use the Hille-Yosida 
theorem to study the higher order regularity of solutions to the "elliptic"equation 
(A — L)w = f, for / G C m ([0, 1]). Finally, in §[[2] we characterize the adjoint 
semi-group, using various abstract results from the theory of semi-groups as well 
as the specific analytic information accrued to this point. This allows us to discuss 
the "forward" Kolmogorov equation, which is crucial for applications in population 
genetics. 

Our focus throughout is on the C° (and C m ) theory, rather than the (simpler) 
L 2 theory. We shall return to a comparison between the C° and I? semi-groups 
elsewhere. 

Remark 1 (Notational Remark). We let I C M be an interval or a ray. We use the 
notation C™ (I) for the space of functions with m continuous, bounded derivatives 
on /; the notation C™(/) for the space of functions with m continuous derivatives, 
and compact support on /. Finally we use C m ([a, b)) to denote the space of func- 
tions with m continuous derivatives on [0, b), tending to zero at b. 
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2 A change of variables 



The first task is to show that by a judicious change of variables, already known to 
Feller [5], any operator of the form © can be reduced to a first order perturbation 
of the exact Wright-Fisher operator Lwf- 

Given the operator ©, assume (by rescaling if necessary) that 



ds 

7T. 







\A0) 



For the heat equation this amounts to rescaling the time variable. We define a 
diffeomorphism of [0, 1] by setting: 

X 

r)(x)\ f ds 



^(x) = sin 2 ^ ) where r/(x 



o 



A short computation shows that L becomes 



where 6(£) has the same properties as the original function b(x), i.e. it is smooth 
on the closed interval, and 6(0) = 6(0), 6(1) = 6(1), so in particular is inward 
pointing at the two boundary points, and vanishes at a boundary point if and only 
if 6 does. 

Henceforth we return to using x as the independent variable and a and 6 as the 
coefficient functions. We assume that our general Wright-Fisher operator has the 
form 

L = x{l-x)dl + b(x)d x , (9) 

where 6(0), —6(1) > 0. As noted earlier, the case of principal interest in popu- 
lation genetics is when these values lie in [0, 1), the so called "weak mutation" 
regime; however, it is no more difficult to treat the general case. In our analysis it 
is actually necessary to consider larger values of these constants, to obtain higher 
order estimate for the cases when |6(x)| lies in [0, 1). 



3 Natural boundary problems 

Because L is degenerate at the boundary of the interval [0, 1], we must carefully 
specify how boundary conditions for L are formulated. The main observation is 
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that operator L := x(l — x)L has regular (Fuchsian) singularities at the two bound- 
ary points x = and x = 1, hence may be studied by standard ODE methods. This 
is something of a red herring, as it is the resolvent operator (A — L)^ 1 that gov- 
erns the behavior of the heat kernel; the resolvents of L and L are fundamentally 
different. The simple vanishing of the coefficient of the second order term in L 
plays a crucial role in the analytic properties of the solutions to the associated heat 
equation and their applications in stochastic processes and population genetics. In 
particular, the diffusion processes associated to L are qualitatively quite different 
from those associated to L. 

We define the notion of indicial roots for L. The complex number s is called 
an indicial root for this operator at x = if 

Lx s = 0(x s ). 

Note that for any s, Lx s = x s_1 c(s, x) where c is smooth up to the boundary, so s 
is an indicial root only if some leading order cancellation takes place. Indeed, 

Lx s = {s(s - 1) + b(0)s)x s - 1 + 0(x s ) 

and this leading coefficient vanishes precisely when 

s{s + (6(0) - 1)) = <=> s = 0,l-6(0). 

There are some subtle differences in the analysis of L when 6(0) < 1, 6(0) = 1 
or 6(0) > 1. The case of principal interest in biology is when < 6(0) < 1 (and 
similarly, > 6(1) > —1). For the boundary conditions we are considering it is 
no more difficult to handle the general case where 6(0) > 0, and, for technical 
reasons, it is actually very useful to do so. The case 6(0) < behaves quite 
differently, mainly due to the loss of the maximum principle. 

It is well known that if 6(0) G R+, 6(0) £ {0, 1,2,.. .}, then any solution of 
Lu = satisfies 

u = m{x) + x 1 ~ m u 2 (x), where m, u 2 G C°°([0, 1)). 

We call the leading coefficients ui(0) and «2(0) the Dirichlet and Neumann data 
of u, respectively. If u solves the equation formally, i.e. in the sense of Taylor 
series, we see that if Uj(0) = 0, then Uj(x) vanishes to infinite order at x = 0; in 
particular, if u 2 (0) = 0, then u(x) G C°°([0, 1)). When 6(0) = 0, the expansion 
has the slightly different form 

u = ui(x) + xlogxu 2 (x), ui,u 2 G C°°([0, 1)), 
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and we again consider u\(0) and 1*2(0) as the Dirichlet and Neumann data, so 
u G C°°([0, 1)) if and only if u 2 (0) = 0. If 6(0) G N, then there is a regular 
solution m, with m(0) non- vanishing, and a singular solution, u 2 of the form 

u 2 (x) =x 1 - b(0) u 2r {x) + log xu 2 i(x), u 2r , u 2 i G C°°([0, 1)), (10) 

with both U2 r (0) and u 2 i(0) non-zero. A solution is regular at zero if and only if it 
is bounded as x —> + . 

There is an almost identical formal analysis of the asymptotics of solutions to 
Lu = f where / G C°°([0, 1]), and for each of these cases we have the 

Proposition 1. Let u be a solution to Lu = f where 6(0) > 0, 6(1) < 0, and 
f G C°°([0, 1]). Then u is smooth up to x = ifu 2 (0) = 0. 

Later in this paper we we also consider the adjoint L f of the operator L, defined 
formally by 

1 f i 

(Lu)vdx = u(L t v)dx, u,v G C°°((0, 1)). 



(11) 







Integration by parts gives 

L*« = ^{x{l-x)u)-l{b(x)u) = l(j- x ({l-x)u)-b(x)u) 
= x(l - x)^u + [2(1 - 2x) - b{x)]4-u - [2 + b'(x)]u. 



For historical reasons, dt — L is called Kolmogorov's backwards equation, while 
dt — L l is called Kolmogorov's forward equation. A short calculation shows that 
the indicial roots of L l at x = are 

s = 0,6(0)-l. (12) 

Suppose that u G C°°([0, 1]). If 6(0), -6(1) £ N U {0}, then we see that 
Jq(Lu)v dx = f Q n(L*u) dx holds if and only if 

hm x _ 0+ [d x (xv(x)) - b(0)v(x)\ = and 
lim^!- [d x ({l - x)v(x)) + b(l)v(x)} =0. 

These limits characterize the adjoint boundary condition. Loosely speaking, a 
function satisfying these conditions is of the form: x fe ^°^ 1 (l — x)^ ( - 1+b ^v(x), 
with v smooth at and 1. In this connection, the boundary conditions for L are 
often formulated as "zero flux" conditions: 

lim x m d x u(x) = 0, lim (1 - x)~ bm d x u(x) = 0. (14) 
We shall return to a more precise discussion of the adjoint operator in § [T2l 
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4 Reduction to model problems on the half-line 



The first step in the construction of the parametrix is to localize the analysis to 
neighborhoods of the boundary points. This leads us to study certain model prob- 
lems, which are the focus of much of the rest of the paper. Indeed, we establish a 
maximum principle, and also derive explicit formulae for the fundamental solutions 
of these model operators, and using these results, prove sharp mapping properties 
for the heat kernels of these model problems. Finally, by perturbation and stan- 
dard parametrix methods, we prove analogous results for general Wright-Fisher 
operators. 

Our goal is the construction of the "heat kernel" or solution operator for the 
initial value problem for the diffusion equation: 

d t u - [x(l - x)d 2 x + b(x)d x ]u = with u(x, 0) = f(x), (15) 

where we recall that an explicit change of variable has reduced the coefficient to 
d x to this normal form. The function 6 is assumed to be smooth on [0, 1], with 6(0) 
and —b(l) non-negative. As above, we let L denote the differential operator 

Lu = [x(l — x)d x + b(x)d x ] u. 

The additional requirements that 

u€C°([0,l] x [0,oo)) nC x ([0,l] x (0,oo))and, 

lim x(l-x)d?u(x,t) = 0fbri>0, (16) 

assure the uniqueness of the solution. This follows from a maximum principle for 
super-solutions. 

Proposition 2. Let u satisfy (fTBT ), and (dt — L)u < in [0, 1] x (0, oo); for any 
T > 0, let D T = [0, 1] x [0, T]. Then 

max u(x,t) = max u(x,0). 
{(x,t)6D T } {*e[o,i]} 

Proof. Let e > 0, then u e = u + e(l + t)~ l is a strict super-solution. The standard 
argument (see [8 ]) then shows that the maximum of u e occurs along the distin- 
guished boundary d' Dt = 8Dt \ [0, 1] x {T}. The regularity hypotheses in (fT6l ) 
show that the maximum cannot occur where x = or 1. For example, if the maxi- 
mum occurred at (0, to), then <9 t u e (0, to) = 0, and d x u e (0, to) < 0. The conditions 
6(0) > 0, and 

lim x(l — x)d x u € (x,to) = 0, 
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contradict the fact that u e is a strict super-solution. Thus 

max u(x,t) < max u e (x,t) = max uJx,0) < max u(x,0) + e. 
{{x,t)£D T } {( x ,t)eD T } {*e[0,l]} {^6[0,l]} 

This estimate holds for any e > 0, which completes the proof of the proposition. 

□ 

The transformation 

y = sin 2 \fx 

maps the interval [0, (tt/2) 2 } bijectively onto [0, 1] , is regular at x = 0, but becomes 
singular at (tt/2) 2 . In the x-variable, the operator Lwf = 2/(1 — y)dy transforms 
to 

Lwf = xd 2 + xc(x)d x . (17) 

Here c(x) is a real analytic function in a neighborhood of x = 0, which includes 
[0, (tt/2) 2 ), with c(0) = -1/3. More generally, L is carried to 

L = xd 2 + xc(x)d x + c(x)d x , (18) 

where c is again a smooth function which satisfies 

c(0) = 6(0). 

The coefficient b of the first order term has the form 

b(x) = bo + xc(x), where bo = 6(0). 

This suggests introducing the half-line model operators 

L\j = xd 2 + bd x . 

In this coordinate the operator L is represented in the form 

L = L bo + xc(x)d x . 

The first order term xc(x)d x is, in a precise sense, a lower order perturbation of 
Lfc . We show, in the succeeding sections, how to construct the solution operator 
for ([T3T ) with regularity conditions (fT6l ), from the solution operators for L^. The 
next several sections are devoted to the analysis of these solution operators. 
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5 Maximum principles and uniqueness for the model op- 
erators 



For the heat equation on Euclidean space, the maximum principle and the resulting 
uniqueness of solutions on R n x [0, oo), satisfying given initial conditions, requires 
that we impose a growth hypothesis on solutions at |a;| — > oo. In this section we 
establish analogous results for the model operators L b , b > 0. 
For any T > 0, we consider solutions on the domains 

D T = [0,oo) x [0,T]. 

Proposition 3. Suppose that b > 0, and 

u G C°(Dt) PI C 1 ([0, oo) x (0,T]) nC 2 ((0,oo) x (0,T]) (19) 

is a super-solution for the heat equation associated to L b , so dfU — L b u < 0. 
Suppose further that 

lim xd 2 x u(x, t)=0,for0<t< T, (20) 

and that for some a, M > 0, \u(x, t)\ < Me ax for all (x, t) G Dt; then 

sup u(x,t) = sup u(x, 0). 
(x,t)eD T xe[o,oo) 

Proof. For A; a non-negative integer, define the function 

V T ,k{x,t) = — ^ ^ k e~ = d^V Tfl . 

Clearly v Tj k(x, t) > 0, is smooth up to x = 0, and one readily checks that (dt — 
Lk)v T ,k = 0. We shall use the v Tj k as barriers. 

The proof is much the same as in the Euclidean case. We first treat the case 
b = k, a non-negative integer. Choose e±, 62 > and r G (0, 1/a), and consider 
the functions 

u e (x,t) = u(x,t) - eiv Tyb (x,t) + — —■ 

These are defined in D T n Dt and satisfy 

8 t u e - L b u e < (21) 

there, hence are strict supersolutions. Our choice of r ensures that v T ^ — ► +00 
more rapidly than e aa; as x — > 00, so u e < for x sufficiently large, uniformly for 
all t G [0, r]; therefore, the standard proof of the maximum principle shows that: 

u e (x,t) < sup u € (x,t) for (x,t) £ D T n D T , (22) 

(x,t)ed'(D T nD T ) 
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where & D T is the distinguished boundary ([0,oo) x {0})U({0} x [0,T)). Clearly 
sup u € is attained at some point in Dt- 

First observe that the maximum of u e (x,t) cannot occur along x = with 
t > 0. Indeed, suppose that the maximum occurs at x = 0, for some to G (0, T). 
If b = 0, then the regularity of u e and (1211 ) give 9 t u e (0, t) < 0, which contradicts 
that u e reaches a maximum at (0, to)- On the other hand, if b > 0, then (|2TI ) would 
imply that 

<9 t ii £ (0, t ) = and 5^(0, i ) > 0, 

which is also a contradiction since it shows that u e (0, to) < u e (x, to) for x small. 
If the maximum were to occur at u(0,T), then dtu e (0,T) > 0, so as before, 
d x u € (0,T) > and hence u e (0,T) < u e (x,T), for small positive x. Hence the 
maximum does not occur along x = 0. Note that this argument relies strongly on 
the assumption of C 1 regularity up to x = 0, as well as ( [20l . 
All of this shows that we can replace (l22l ) by 

u e (x,t) < sup u e (x,0) < sup n(x,0) + €2. 
xe[o,oo) xe[o,oo) 

Now letting e±, £2 — ► 0, we conclude that 

sup u(x,t) < sup u(x,0). (23) 

{x,t)&D T r\D T xe[o,oo) 

If t > T, then this completes the proof of the proposition; otherwise repeat the 
argument recursively, with u\{x, t) = u{x, t + t) replacing u(x, t). Finitely many 
iterates produce the desired conclusion. 

We now turn to non-integral values of b. If k < b < k + 1, for a non-negative 
integer k, then 

(d t - L b )(u - eiv T , k +i + Y~jT^) = 

(d t - L b )u + e X (b- (k + l)K,fc+2 " 7T~7\2 < °- 

The inequality uses that tv,fc+2 > 0. We can argue exactly as before with 

u e (x, t) = u- eiv T)k+ i + , 

to obtain (l23l . and iterate, if needed, to complete the proof of the proposition. □ 
The uniqueness of classical, "tempered" solutions is an immediate corollary: 
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Corollary 1. Let u satisfy d 19b , (1201 ), and be a solution to dtu — L^u = Ofor some 
constant b > 0. Ifu(x, 0) = and 

\u(x,t)\ < Me ax for (x,t) £ D T (24) 

for some a, M > 0, then u = in Dt- 

Remark 2. The discussion of indicial roots in Section [3] shows that other solutions 
to the initial value problem for dt — Lb, have an asymptotic expansion at x = 
with the term x 1-6 , so even if the initial data is smooth, those solutions are not C 1 
up to x = and do not satisfy (l20l . Following Feller, the local boundary condition 
that singles out the smooth solution is the zero flux condition: 

lim x b d x u(x, t) = 0. (25) 

6 Fundamental solutions for the model operator 

Consider the heat equation 

dtu = LbU, x>0, u(0,x) = f(x) (26) 

with boundary condition at x = for t > dictated by the demand that u(t, x) G 
C m up to x = for all m G N provided the initial condition / also lies in C m (and 
has moderate growth as x — > oo). This is clearly needed for the solution operator 
to (l26l > to define a semi-group on C m . 

The goal in this section is to find an explicit expression for the fundamental 
solution of this problem. Our strategy is as follows: We first study the Fourier 
representation of the heat kernel for Lq\ this turns out to be somewhat simpler to 
analyze, and it suggests the general form of the corresponding kernel when b > 0. 
It is useful to have both a Fourier representation and the explicit Schwartz kernels 
of these operators, so we include material describing both approaches. 

Remark 3. We would again like to thank Charlie Fefferman for showing us his 
construction of the kernel k^(x,y), essentially (|28T ). on which all our subsequent 
development is based. This kernel also appears in O. 

6.1 Fundamental solution for d t — L 

We begin by seeking a function E (x, £, t), which is a solution to dtE° = LqE°, 
with initial condition E°(x, £, 0) = e lx ^. A first guess is that E° should involve a 
Gaussian, but because of the natural appearance of the variable ^fx in this problem 
(or, equivalently, the fact that the dilation (t, x) h-> (At, Ax) preserves the equation 
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under consideration), we are led to the ansatz E°(x, £, t) = exp(x0(^, t)) for some 
as yet unknown function (p. Inserting this into the equation leads to the initial value 
problem 

d t <p = <p 2 , <H£,o)=^£, 

which yields that 

As an oscillatory integral, the corresponding Schwartz kernel is 

1 f'°° 

k°(x,y) = —J £°(x,e,t)e-^de. (28) 

With a bit of algebra, one sees that the integrand decreases exponentially as a func- 
tion of Im^ when x,t > and y < 0. Thus, a contour deformation shows that 
k® vanishes in this region. In other words, k®(x, •) is supported in y > when 

x,t > 0. 

For the special case 6 = 0, the condition u(0, t) = defines a CQ-semi-group, 
which turns out to be a sub-semi-group of that defined by kf. We call this the 
"Dirichlet semi-group." Since fc°(0,y) = 5(y), the kernel we have constructed is 
not the Dirichlet heat kernel, but can be obtained by a simple modification: 

k° t ' D (x,y) = k° t (x,y) - e~ x N(y). (29) 



In other words, the fundamental solution for (1261) with Dirichlet boundary condi- 
tions at x = is 

A^( X) y) = _L y°° ( e («*)/(i-«0-itf _ e -*/*-itf) ^. ( 30) 

The integrand here again diverges like l/£ as |£| — > oo, so we interpret (l30l) as 
an oscillatory integral via the regularization (obtained by a formal integration by 
parts) 

oo 

*< (X ' y) = -Uji (tC + *) 2 2^ ^ 



This satisfies the same heat equation, but now vanishes x = 0, so is indeed the 
Dirichlet heat kernel. 

It is possible to find an explicit expression for k®' D in terms of elementary 
transcendental functions. First note that k t ' is a function of 

x a a y 

a = — and p = — . 

t 1 t 
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Indeed, setting rj = ££, and z = rj + i, we find that 

fcO.D, ^.M f ev+ie-^ dr, _ {a+(3) ( a\ f e~~^ dz 
kt {X ' y) ~ W) J (r? + .) 2 2vr" W)J z 2 2vr' 

-oo Ti 

where, by definition, T T = {z : Imz = r}. By an elementary contour deforma- 
tion, 

f e'f-' 1 ^ dz _ r e^-W* dz r e^^ 2 dz 

J z 2 J z 2 J z 2 2^ 

ri r T \ z \=i 

for any r < 0. The first term on the right vanishes since it can be made arbitrarily 
small by letting r — > — oo, so it suffices to compute the integral on the unit circle. 
For this, note that the coefficient of 1/z in 

oo ^ oo ^ 

z~ 2 exp(ia/z) x exp(—0z) = z~ 2 ^ —{ia/z) j ^ — {—i(3z) k 
is equal to 

oo -, oo / a\£ 

where 1\ is the modified Bessel function of order 1. This yields, finally, the explicit 
formula 

u o,D, y. 1 -Eta fx {2y/xy\ 

for the Dirichlet fundamental solution of the model heat equation. It is useful to 
represent this kernel in an alternate form. For b > we define the entire functions 



3=0 

An elementary calculation shows that ^ satisfies the ordinary differential equation: 

+ bip' b - ^ b = (33) 
The Dirichlet heat kernel is then given by 
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The kernel k® can be expressed as 

k° t {x,y) = e-T5(y) + k^ D (x,y). (34) 

The relationship between tp2 and the /-Bessel function implies the asymptotic 
expansion: 



47TZ4 



oo 



See H 8.451.5] 



6.2 Fundamental solution for d t — L b 

We now undertake a similar analysis of the fundamental solution for the problem 
(l26l ) for < b. As explained in §[3]the boundary condition at x = which should 
guarantee that solutions are smooth up to x = is the analogue of the Neumann 
condition, i.e. the one which excludes the term x 1 ^. We denote by k b {x,y) the 
heat kernel for this problem with this zero-flux conditions. 

As before, we first determine the Fourier representation of k\ using the ansatz 
that E b (t, x, £) = ^(t£)e -x ^*'^; this leads fairly directly to the expression 

E b (t,x,0 = (l-itO- b e^, (35) 



and hence, (when b < 1) as an oscillatory integral, 



* S(l - v) -s/ < »^W* <36) 

If / € C°([0, oo)) has an absolutely integrable Fourier transform, then this formula 
can be interpreted to mean that 

OO OO ixj ^ 

k^, v)mi y = lJ^J^. (37, 

-oo 

In order to evaluate (l36l ). we start off as before, setting a = x/t, (3 = y/t, and 
expressing (l36l > as a contour integral over Im z = 1 : 

nib R 2+i 

-Rl+i 
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To define z~ b , we take aigz = on the positive real axis and consider z~ b as 
defined on the plane cut along the negative imaginary axis. 

Now change variables, setting z = \Ja/ [5t, and for simplicity write £ = 
y/aP, to write 



A^ z )z- b dz 



Imz=l V @ Jlmr=yJ (3/a 

Since the integrand decays as | Rer| —>■ oo and as Imr — ► — oo, we can deform 
the contour to yield 

r 37T — 7T, 



A = {ye— : y £ (oo,l]}U{e 1, 



€ ly,— ]}U{ye 2 :y e [l,oo)}. 



Inserting the explicit parameterizations of the various parts of this contour into the 
integrand e^(~~ T ) % and simplifying yields the absolutely convergent representa- 
tion: 



7T OO 

J e 2Ccos *cos(6- l)0#-sinvr(6- 1) J e 



,b ' 



Changing the variable to y = e* in the second integral on the right, we get 



-2Ccoshi e -t(6-l)^_ 



2 y e 2 ? cos ^cos(6- l)0#-2sin7r(6- 1) ^ e~ 



This combination of integrals appears as Formula 8.431.5 in |6j, and is seen to 
equal 27r/b_ 1 (2£). Putting all of these calculations together gives the two equiva- 
lent expressions 



and 



t \y 
k° t {x,y)dy = 



(38) 



(39) 



b -z+y , (xy\ dy 

tj 1 Ht) 7 

Using the classical asymptotics for we see that as z — > 00, ipb has the asymp 
totic expansion 



1 + E £M 
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(40) 
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From these explicit formulae it is evident that the kernels k\{x,y) are pointwise 
positive. Once it is shown that they generate semi-groups, it follows that these 
semi-groups are positivity improving. The Fourier representation has a very useful 
consequence: 

Lemma 1. /// G Cf([0, oo)) satisfies dif(0) = OJor 1 < j < m, then 

oo oo 

di J k h t [x, y)f(y) dy = J k\ +3 {x, y)cP y f(y) dy, for l<j<m. 
o o 



Proof. This follows easily from (I37l >. as the hypothesis implies that we can dif- 
ferentiate under the integral sign. The vanishing conditions on dif(0) give that 
£™~l/(£) e so that 

%f(0 = (iO j f(0, forl<j<m-l. 

If £ m / (£) £ -L 1 (M) , then we can immediately extend the formula to j = m. In 
general, approximate / in the C m -norm by a sequence {/„} C C m , which satisfies 
£>if(0) = 0, for 1 < j < m and this additional integrability hypothesis. A 
limiting argument then establishes the truth of the lemma for j = m in the stated 
generality. □ 

Remark 4. If < b < 1, then, as shown in 0, one can define the Dirichlet 
problem for <9 t — L^. For this range of parameters this corresponds to having a non- 
zero coefficient of x 1-6 , hence a non-smooth solution. For b > 1, Feller showed 
that (|25T ) is the unique local boundary condition that defines a positivity preserving 
semigroup. 

6.3 Representation formula 

As a first application of the existence of these fundamental solutions, we prove a 
representation formula for the solution to the initial value problem 

dtu = Li,u with u(x, 0) = f(x), (41) 

where / and u satisfy ((24)) and in addition, 

lim x b d x u(x, t) = for b > 0, lim xd x u(x, t) = if b = 0. (42) 

Note that if u is in C 1 {Dt), then it automatically satisfies this boundary condition. 
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Proposition 4. Suppose that u £ C°(D T )nC 2 ((0, oo) x (0, T)) satisfies (|4TI (02]). 
7/"/, u, dtu and d x u satisfy M4\ then for all t < a , u is given by the absolutely 
convergent integral: 



i(x,t) = J k b t (x,y)f(y)dy. 



Proof. For b > 0, the proof of this proposition is a simple integration by parts 
argument. It uses the fact that k\ also satisfies the adjoint equation and boundary 
condition: if t > 0, then 

9 t kt{x, y) = (ply - bd y )k\{x, y) and lim [d y yk\(x, y) - bk b t {x, y)\ = 0. 

y-*o+ 

Moreover, it is not difficult to show that as y — > oo, with x in a compact set, we 
have 

26-3 y 

(43) 



k$(x,y) <Cy * e t 
With these preliminaries, we can integrate the equation satisfied by u to obtain 

t — € OO 

= J J k h t _ s (x, y) (d s u(y, s) - L b u(y, s)) dyds. 
o o 

Here e is a small positive number. Provided t < a -1 , the estimates satisfied by 
u, dtu, d x u, and (l43l justify the subsequent manipulations of this integral. 
We first integrate by parts in s to obtain: 



[u(t - e, y)k b e (x, y) - u(0, y)k b t (x, y)\ dy 



oo t— e 



L b u(s,y)k b _ s (x,y) - u(s,y)d t k b _ s (x,y) 



dsdy. 



o o 



Now integrate by parts in y using the boundary conditions satisfied by u and k\_ 
at y = 0, and the estimates (l24l and (l43l . to conclude that since t < a^ 1 , then 

oo 

J [u(t - e, y)k b (x, y) - u(0, y)k b (x, y)]dy = 



oo t— e 



u(s, y) L^ t fe _ s (x, y) - d t k b t _ s {x, y) 



dsdy. 



o o 
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For any e > 0, the right hand side of this equation vanishes identically, and there- 
fore 

oo 

*-.,,«(,,,)*- MM«(..v>*. 



Letting e — » gives the result. 

For b = we proceed a little differently If we write u = Uq + /(0), then uo 
also satisfies (e?t — Lq)uq = 0, Uq(0, t) = 0, and the same regularity conditions as 
those satisfied by u. Arguing as above, 



u (x,t) = I k°> D (x,y)(f(y)-f(0))dy, 



and hence 

u(x,t) = I kt(x,y)f(y)dy. 



□ 

As a special case we can demonstrate that the kernels {k\ : t > 0} have the 
semi-group property. 

Corollary 2. Ift, s are positive numbers and b > 0, then 

k t+s( x ,v) = J kt(x,z)k b s (z,y)dz. 


Proof. For fixed y, the kernel k b (-,y) satisfies (9 t fcj(x,y) — Lf,k^(x,y) = 0, and 
is C°° in [0, oo) x (0, oo). This solution decays exponentially as x — > oo, so the 
maximum principle implies s i— > k^ +s (x,y) is the unique smooth solution of this 
PDE with respect to the variables (s, x), with initial data k\ (-, y). The proof now 
follows from Proposition (0]). □ 

Remark 5. A similar argument using the uniqueness of the solution to the Dirichlet 
problem when 6 = leads to a proof that 

oo 

0,D, n _ f U 0,D, \ ,.0,1) 



k»£(x,y)= k^(x,z)k^(z,y)dz. 
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It is also true that if / € C£°([0, oo)), then 

oo oo 

I hk a (x,v)f(y)dy = J k° t (x,z) 



We conclude this section with: 



k° s (z,y)f(y)dy 



dz. 



Proposition 5. If f G C ([0, oo)) /za^ arc absolutely integrable Fourier transform, 
then 



lim 

6^0 



k b t {x,y)f{y) dy 



kt(x,y)f{y) dy. 



This convergence is obvious using the Fourier representations d28l ) and 051 ) of 
the kernels and fcj! 1 . If / is absolutely integrable, then, for any fixed < T the 
convergence is uniform on [0, oo) x [0, T\. 



7 Mapping properties for the model operators 

We can now prove the key fact that the diffusions associated to the operators 
are semi-groups, which, for every m £ N U {0}, preserve the spaces C m ([0, oo)), 
see Remark [6] We prove this in two steps. The first is to show that the kernels map 
polynomials to polynomials; then, decomposing an arbitrary C m initial condition 
into a polynomial (its Taylor series) and a remainder term that vanishes to order m, 
the fact that the solution with initial condition given by this remainder term is also 
C m follows from the maximum principle. 

Suppose that the initial data is a polynomial, f(x) = X]j=o a j x ^- It i s clear 
that the formal exponential 

e=o l - 

makes sense and in fact is a polynomial which is the unique moderate growth 
solution of (|4TT) . We simply observe that L^x J is a constant multiple of x^~ e , hence 
this vanishes as soon as I > j. Thus the sum on the right above is finite and is a 
polynomial in (x,t). It solves the equation by the usual elementary calculation. 
Therefore, by the uniqueness theorem, this is the only (exponentially bounded) 
solution of this problem, and hence, by the representation formula, we must also 
have that 

o 

We now turn to the case of general C m data. 
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Proposition 6. Suppose that f £ C m ([0, oo)) has compact support. Then 

poo 

u(x,t)= / k b t (x,y)f(y)dyeC°([0,oo) t ;C m ([Q,oo) x )), 
J o 

solves the initial value problem (|26l >. The norm of the difference, \\u(-, t) — f\\c m 
tends to zero as t — > 0, and therefore k\ extends to define a C m -semi-group for 
each m = 0, 1, 2, 

Remark 6. We let C m ([0, oo)) denote the closed subspace of C m ([0, oo)) consisting 
of functions / with 

lim cP x f{x) = for < j < m. 

x— »oo 

This is the closure of C™([0, oo)) in the C m -norm. 

To prove the proposition, the following two lemmas are useful 

Lemma 2. For b > and < i < L, there is a constant Cb so that if f € 
C£([0, L)), then, for x > L, and t < 1, we have the estimate 

oo 

j k b t (x,y)f(y)dy 
o 

If f £ C°(/,oo), is of tempered growth, then for x < I, andt < 1, 

e t 

\u(x,t)\ <C btN \\f\\( N) - b . 

Here N € N is chosen so that 

\\f\\ {N) :=sup(l + \x\)- N \f(x)\<oo. 

0<x 

The proofs are elementary using the asymptotic expansion for ipj, and are left 
to the reader. It is also useful to have the following 

Lemma 3. Ifb > 0, and f £ C™((0, oo)), then 

oo 

u(x,t) = J k h t {x,y)f{y)dy, 
o 

solves d26l ) and 

lim \\u(-,t) - f\\ C m = 0. 



< c s ||/ll= 



L\ L - 



4 e" 
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Proof. It is clear that dtu — L^u = in [0, oo) x (0, do). Suppose that supp / C 
[l,L]. The previous lemma shows that for any r\ > 0, the functions u(-,t)(l — 
X[l-r),L+rj\(')) converge uniformly to zero as t — > + . If r\ < 1/2, then xy > l 2 /2, 
for (x,y) G [I — rj,L + rj] x supp/, and we can therefore use the asymptotic 
expansion for tjjj, to conclude that t)x[i-r),L+rj\ (') converges uniformly to /. As 
/ G C™((0, oo)), we know from Lemma Q] that 

oc 

d 3 x u{x,t) = j k b t +J (x,y)dU{y)dy, for j < m. 



o 



The argument above applies to show the uniform convergence of the derivatives 
diu(-,t) to dif, for 1 < j < m. □ 



Proof of Proposition [6] Fix m and / £ C™ [0, oo) . From the expressions (l40l) and ((34 
for fc£ it is clear that, if / e C°, then u G C°° ([0, oo) x (0, oo)), and d t u - L b u = 
in [0, oo) x (0, oo). Choose a smooth cutoff function x( x ) which equals 1 near 
x = 0, \x( x )\ — 1. an d which vanishes for x > 1. Let 

^ /0-)(Q) 
i=o J ' 

be the Taylor polynomial of order m for /; thus 

/(») = /(») - X(x)q(x) = o(x m ), asx\0. 
We observe that 

/•oo 

u xq (x,t) = I k b t (x,y)x{y)q{y)dy 

OO rOC 

kt(x,y)q(y)dy - / k\ (x, y)(l - x(y))q(y) dy. 
io Jo 

By the remarks above, the first term on the right is a polynomial, and therefore in 
C°°([0, oo) x [0, oo)). On the other hand, (1 — x{y))q{y) i s supported away from 
y = 0, and it is elementary from the regularity properties of k\ proved in the last 
section and Lemma |2] that 



o 



oo 

b 



k b t (x,y)(l - X (y))q(y) dy £ C°°([0,oo) x x [0,oo) t ), 



o 



and tends locally uniformly to (1 — x( x ))q{ x ) as t ^ + . Together with LemmaO 
this shows that u xq (x, t) tends to xq i n C([0, oo)) for all I G N. Thus it remains to 
prove that 



POO 

u(x,t) = / kt(x,y)f(y)dy 
Jo 
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is also in C m and tends to / in the C m -norm as t — > + . 

First consider the case m = 0. For every e > we can choose 5 > so that 
f(x) < e when x < 5. Decompose f(x) = x{x/S)f(x) + (1 — x( x /S))f(x). 
Applying k\ to the second term gives a smooth function v(x,t) in C°°([0,oo) x 
(0, oo)), which satisfies d t v — L^v = 0. Lemma [5] shows that v G C°([0, oo) x 
[0,oo)), and tends uniformly to (1 — x( x /$))f as * — * + . As A;* is pointwise 
positive, and has integral 1 for all t > 0, we see that 



oo 

6/ 







< e. 



K{x,y)x{y/5)f{y)dy 

Since e is arbitrary, this shows that 

limsup ||u(-,t) - f\\c° = 1. 

If m > 0, then dif(0) = 0, for 1 < j < m, and we can apply Lemma[T]to 
conclude that 

oo 

d J x u(x,t) = J kt +3 (x,y)d 3 y f(y)dy, for 1 < j < m. 
o 

The C° -argument then applies to complete the proof. □ 

As a corollary of these results we can extend Lemma Q] to general data in 

C c m ([0,oo)). 



Corollary 3. Ifb > 0, and f G C™([0, oo)), and 





oo 

^u(x, t) = J k b t +j (x, y)d> v f(y), for 1 < j < m. (44) 
o 

Proof. The Proposition implies that u is a solution to dtu — L^u = 0, which 
belongs to C m ([0, oo) x [0, oo)), and satisfies 

lim u(-,t) = f, 



22 



with convergence in the C m -norm. One can therefore differentiate the equation 
satisfied by u to conclude that 

dtdlu — Lb+jU = for 1 < j < m, and lirri^ &> x u = d^f. 

For each j, the right hand side of (l44l ) is another solution, uj to this initial value 
problem, that also satisfies the hypotheses of the maximum principle. Thus Uj = 

diu for 1 < j < m. □ 

We can use this result to study the regularity properties of the operator 

t oo 


A case of particular importance in applications is when g G C°°([0, oo) x (0, oo)), 
for example g is any solution to (l26l ). We begin with a lemma. 

Lemma 4. Suppose that b > 0, and f G C™([0, oo)); then 

oo 

.(„«)=/*(,„)/(,)* 



oo 

d{u(x,t) = j k b t (x,y)Llf(y)dy 
o 

/or 2j < m t > 0. 

Proof. Suppose that m = 2 and / vanishes to order 2 at j; = 0. Using the fact that 
(dt — (LyY)kt(x, y) = and a simple integration by parts argument we easily see 
that 

oo 

d t u(x,t)= / k b t {x,y)L b f(y)dy. 



o 



The second order vanishing of / at implies that the boundary terms at x = 
vanish. 

Choose a smooth function x, which equals 1 in [0, 1] and is supported in [0, 2]. 
If / does not vanish at x = 0, then we let 

<fc(s) = X(x) (/(0) + /'(0)* + /"(0)y) , and / = / - g 2 . 



23 



From the argument above, 

oo 

d t J k b (x,y)f(y)dy = J k b (x,y)L b f(y)dy, 
o o 

and then by Proposition [6l 

oo oo 

d t J k b (x,y)q 2 (y)dy = J k b (x,y)L b q 2 {y)dy. 
o o 

This argument can be applied inductively to obtain the case of general m. □ 

Combining this lemma with Corollary [3j we conclude 
Corollary 4. Suppose that b > 0, and f G C™([0, oo)), then 



u(x,t) = I k t (x,y)f(y)dy, 

satisfies 



o 



%d*u(x,t) = J k b t +k (x,y)L{ +k d k y f(y)dy, 
o 

provided 2j + k < m and t > 0. 



We can now examine the regularity of K b g. 

Proposition 7. Ifb>0,T>0, and 

g G C 6 °([0, oo) x [0, T)) D C b m ([0, oo) x (0, T)), 

/or < t < T, the derivatives cP t d\K\g are continuous for 2j + k < m; 
moreover, ifm > 2, then 

(d t - L b )K b g = g (45) 
Proof. For t > 2e > 0, we define 

t—e oo 

K t,e9( x ) = J J k b _ s (x,y)g(y,s)dyds. 
o o 
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For e > it follows from the corollary, that the derivatives cf t d^.K\ t g exist, pro- 
vided 2j + k < m and can be expressed as 



i 

2 OO 



4d k x Kl e g{x) = J Jdld k x kl s (x,y)g(y,s)dyd S + 

k b +*(x, y)L{ +k d k y g(y, s) dyds. 



o o 

t — £ OO 



i 

2 



If 2j + k < m, then & t d%K$ e g(x) converges locally uniformly to a continuous 
function. This establishes the existence of these derivatives. If m > 2, then differ- 
entiating, for any < e < t we have that 

OO 

(d t - L h )Kl S = J k b e (x, y)g(y, t - e)dy. (46) 
o 

As m > 2, the limiting function K\g has one time, and two spatial derivatives, 
which are the limits of the corresponding derivatives of K b e g. Thus letting e — > + 
in (06]> gives (|45l) . □ 



We can finally extend the convergence result Proposition [5] to arbitrary data in 

C°([0,oo)). 

Proposition 8. Fix f E C°([0, oo)), then 



roo poo 

Jim/ kt(x,y)f(y)dy = k^(x,y)f(y)dy. 
Jo j 

For any T > convergence takes place in the C°([0, oo) x [0, T])-topology. 

Proof. Choose a sequence {/ n } C C)?([0, oo)) with all /„ G L 1 , and such that /„ 
converges uniformly to /. Let u b denote the solutions to (l26l ) with initial data / 
and u b n the solutions with initial values f n . The maximum principle implies that, 
for any n, 

II"" 6 _ "°llc°([0,oo)x[0,T]) < 2 II/ ~ /nlloo + \\u n - Unllc°([0,oo)x[0,T])- 

Given e > 0, fix some n so that ||/ — / n ||oo < £• Applying Proposition [5] to the u b n 
gives a positive bo so that if 6 < bo, then 

ll u n - u n\\c°([0,oo)x[0,T}) < € , 

as e > is arbitrary, this completes the proof of the proposition. □ 
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8 Perturbation estimates 

In order to use the model heat kernels k\ in perturbative constructions for the heat 
kernels of general Wright-Fisher operators, it is necessary to prove estimates in C e 
for every £ > for operators of the form 

t oo 

g i-> A\g{x) = J J kf_ s (x, z)h(z)zd z g(z, s) dzds, (47) 
o o 

were h G Q?°([0, oo)) is fixed. 

Corollary [3j the Leibniz formula and the mapping results established in the 
previous section show that if g G C e , then 

diA b t g = A b +%g + 

£ t oo 

SO)/ Jk^(x,z)di(zh(z))di +1 -ig(z,s)dzds. (48) 
i =1 o o 

Hence to estimate A\g on for arbitrary I > it suffices to prove mapping 
properties for A\ on C° for all b' > 0. 

The following lemma is the key to all that follows: 

Lemma 5. There is a constant Cb, defined for each b > and uniformly bounded 
for b in any compact interval [0, B], such that if f G C?([0, oo)) and 



u(x,t) = J k t (x,y)f(y)dy. 
o 

then 

\d z u(z,s)\<C b M^. 

S + y/ZS 

Proof. It is enough to prove that for each b > there is a constant Cb so that 



<t>s,b{z) :-- 



b 







S + x/SZ S (l + y/z/s) 



and that Cb is uniformly bounded above on any interval (0, B], i.e. it does not 
depend on a positive lower bound for b. The case b = is then obtained from a 
separate limiting argument. 
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-) < (pr) - 1* (pr 



We first compute that 

00 



Set w = y/s and A = z/s, so that </> 8t b(z) = (f ) , where 

00 

<p(X) = J w b - l e- w e- x \wij' b (\w) - ^ b {Xw)\dw. 


Hence we only need prove that tp(X) < C b /(1 + V~X). 
Using the asymptotic formulas 

y;4~2e 2v/ "' _i 

ipb(w) ~ p= — (l + 0(w 2)) and 

/47T 



2/ ' 



-(i+ileV 5 1 
-^(l + 0(.-.)), 



we see that <p(A) < C^a f° r A < A. It is less obvious that this constant is bounded 
as b — > since the integrand appears to become nonintegrable at u; = in that 
limit. To analyze this, define ipb( z ) = 4>b{ z ) — V'b(O) and recall that Vfe(O) = 
l/r(6), so that 



<p(A)< 



..6-1. 



|w^(Aw;)| + |^ 6 (Au/)| + 



1 



r(6)J 



dw. 



The integral of the last term in brackets is identically equal to e _A . As for the 
other two terms, note that the coefficients in the error terms in these asymptotic 
formula? are bounded as b — > 0, so the integral from 1 to 00 converges uniformly, 
independently of b. If w < 1 (and A bounded), these terms are 0(w), hence this 
part of the integral is also uniformly bounded. Hence C bj \ is uniform in b G [0, B] 
for any fixed B, A. 

Now consider what happens when A — > 00. Suppose first that b > 0. Break the 
integral defining ip into the sum J' b + J' b ', where J' b is the integral from to 1/ y/X 
and J b is the integral from 1/y/Xto 00. It is straightforward that J b < Ce~ cX , for 
c, C > which are independent of 6. For the other part use the asymptotics of ip b 
and if/ b to get 



dw 
w 
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Changing variables to y = yw/X — 1 transforms this to 

OO 

J' b '<C b V\ j (l + y) b -h- x y 2 \y\dy, 



-i 



which now, by Laplace's method, satisfies J{j < -^j=. This proves the estimate for 
b > bounded away from oo. 

To finish the proof, observe that for any 7] > and fixed z, s > 0, we have 



lim 

6^0+ 

V 



oo oo 

J \d z k b s (z,y)\dy = J \d z k° s > D (z,y)\dy 



Since the constant Cb is uniformly bounded as b — > 0, there is a constant Co so that 
for any rj > 0, 



>3C 



u.D, ,i , , Co 



IWr(«,y)l<ii/< 



'/ 



8(1 + V^/s) 



The right hand side is independent of r\, so letting t] —> gives the same estimate 
for the integral on all of M + . Setting z = shows that 



| \{d z k° s > D )(Q,y)\dy 



^2(0) 



Since 



fc s °(z,y)/(y)dy = y k° s > D (z,y)(f(y) - f(0)) dy + f(0), 
o o 

this completes the case when b = as well. □ 
Mapping properties of (A b t ) j 

Now we turn to estimates of A h t and its iterates. Because the proofs of the next two 
Propositions are quite technical, we state the results here and relegate their proofs 
to Appendix |A] at the end of the paper. 

We assume that g G C°([0,oo) x [0,T]) n ^((OjOo) x (0,T]) with a very 
specific blowup for d x g as x — > + , as suggested by Lemma|5] 
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Proposition 9. Define the sequence of constants 



( \2±i 



r 



3+1 
2 



For any smoof/z h with supp /i C [0, L], if g G C°([0, oo) x [0, T]) n C 1 ((0, oo) x 
(0, T]) satisfies 



f/ien 

1(4)^(^)1 < ^MCr'CVZ^IUyti (50) 

and 



w/iere Cf, « constant appearing in Lemma\5\ 
Remark 7. The hypotheses of this proposition imply that 

\yh{y)dy9{y,s)\ < - — 7= — , 

V s 

so we could replace k®_ s by k^!° s in the definition of A®. With this choice, A®g(0) = 
0. However, k^' D does not satisfy Corollary |4j and its use would also complicate 
the derivation of the estimate for diA® g{x), so it is simpler to use the kernel k®_ s 
in the definition of A® . 

For the higher norm estimates for A\ and its iterates, it is useful to introduce, 
for T > and I G N, the norms 

IMIc*>°°[0,T] = ™^,ll^(-»*)llc<([0,oo)) 



The maximum principle and (1481 ) immediately give the 

Lemma 6.1fge C l h ([0, oo) x [0, T]), and 1 < p < £ - 1, j G N, f/ierc 
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Proposition 10. Define the constants Dj inductively by 

and Dq = 1, so that Dj < C f +i for all j > 0. Let h be any smooth function 
with supp h C [0, L\ and suppose that the function 

<?eC £ ([0,oo) x [0,T])nC m ((0,oo) x [0,T]), 

satisfies 

\di +1 g(x,t)\<^= (52) 
x/xt 



for some constant M > (depending on g). Then there are constants CT,L,e,b, C' T L ib 
so that, for j £ N, we have 

\di[A\Y g{x)\ < — ^(M + |bllc^[o,r])(^,L,^ll^ll^) J ^, (53) 

and 

< ^( M+ ''^'^[°.^^'' fe ''^) J '^. (54 ) 

Off-diagonal estimates for (A b t ) j 

The final set of estimates we derive for the operators A\ involve the off-diagonal 
behavior of the Schwartz kernel of the infinite sum 

00 00 rt nsi roo poo 

Y J (A b t yk b t (x,y) = Y, / ••• / / •••/ A t- Sj (x,zi)x 

A^.^.^iz!, z 2 )... A b S2 _ Si fa-!, zj)k b Sl ( Zj ,y) d Sl . . . d Sj d Zl . . . d Zj . 

The most precise bounds are best described in a blown-up space. For our purposes, 
the slightly cruder estimates derived here suffice. 

We begin by observing that if \x — y\ > a > 0, then k\ (x, y) < Cy b ~ 1 e~ c l t 
(where c is any number less than (y/x — yjy) 2 )- In fact, k b (x, y) = y b ~ 1 F(t, x, y) 
where F is smooth away from the diagonal x = y but up to x = and y = 0, and 
satisfies |F| < Ce _C// * where c depends only on \x — y\. This is proved using the 
explicit expression for k\ , and separately considering the behavior in the regions 
xy > t 2 and xy < t 2 . Using this we can now prove the 
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Proposition 11. Fix a > 0, and k, £ € N; then there exist constants C,c,B > 
depending on a, \ \h\ \ 00 , L, k and i such that if \x — y\ > a and < x, y < L, then 

\d k x {yd y f[{A b t yk b t ){x,y)\ < CD 5 Bie-^y h -\ 

where Dj are the constants appearing in Proposition [70] (These may be replaced 
by the constants djfrom Proposition^when k = £ = 0.) 

Proof. As just indicated, this assertion is clear from the formula for k\ when j = 0, 
so we proceed by induction, assuming that we have proved it for all powers up to 
some j. 

First use a smooth partition of unity to decompose k\ = k\' h + k\ ' b where k\ ,h 
is supported in the region \x — y\ < a/2 and y 1 ^ b k^ is C°° when x, y, t > and 
satisfies |y 1 ~ 6 A;^'''(x, y)\ < Ce~ c l l for all x, y. Using Proposition |9](or Proposition 
[l0]for the higher derivatives), we have that 

\(A b y +1 k^ b (x,y)\ < CdjB^e-^y*- 1 . 

To estimate the other term, let us begin by writing 

A\ o k]' b (x,y) = / / k b _ s (x,z)h(z)zd z kl' b (z,y)dzds 



o Jo 

t/2 roo 

/ -(zd z + l)(kl s (x,z)h(z))kl> b (z,y)dzds 
'o Jo 

ft fOO 

+ / k b _ s (x,z)h(z)zd z kl' b (z,y)dzds. 

Jt/2 JO 



The integration by parts used to obtain the first term on the right is valid because 
zk b _ s (x, z)kl ,b (z, y) < Cz b and h(z) has compact support. For this first term, use 
that \x — z\ > a/2 and \z — y\ < a/2 to get 

\{zd z + l)k b s (x,z)\ <Ce- c ^z b -\ and \kl' b (z,y)\ < Cs~ 1 / 2 y b ~ 1 , 

which shows that this term is bounded by B\e~ c l t y b ~ 1 . The second term is bounded 
similarly, using 

\d z k\\z,y)\<B 2 s-^ 2 y b -\ 
Taken together, we obtain that 

lAlk^y^KBe-^y"- 1 

where B depends only on the quantities indicated. Finally, use Propositions [9] and 
[JO] to obtain 

\{A\y o{A\kl> b )\<Cd j Bi+ x e- c l t y h -\ 
as claimed. The estimates for the higher derivatives are proved similarly. □ 
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From these estimates we now obtain the 
Corollary 5. The kernel q^(x, y) of the operator 



oo „ 

/~$/ = £(4) j J k b t (;y)f(y)dy 



3=0 

can be written as 

$ = y h - l q^ Cg (x,y), where ^' rcg G C°°([0, 1] X [0,1] x (0,oo)). 

For each d > 0, there is a constant > 0, so that in the off-diagonal region 
{(x, y) : \x — y\ > d}, we have the estimate 

\q t \x,y)\<Ce-¥y b - 1 , (55) 
with analogous estimates for the derivatives d^(yd y ) q t (x, y). 

9 Construction of the heat kernel for a general Wright- 
Fisher operator 

After this long excursion into the analysis of the model problems we are now pre- 
pared to construct the heat kernel for the full generalized Wright-Fisher operator 

L = y(l-y)d 2 y +b(y)d y . 

Let 

y = sin 2 yfxi and 1 — y = sin 2 ^/x^. 

According to the discussion in Section [4] pulling back L to these coordinate charts 
gives 

L t = x l d 2 Xi + b d Xe + x e ci(x e )d Xe and L r = x r d 2 Xr + hd Xr + x r c r (x r )d Xr , 
where 



6 = 6(0) and b x = -6(1). 

that u solves { 

the functions 



Suppose that u solves ( fT5T >, with u(y, 0) = f(y). Then on the interval [0, (f ) 2 ) 



ug(x£, t) = u(sin 2 ^/xi, t) and u r (x r ,t) = u(l — sin 2 \fx~ Tl t) 
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satisfy 

(d t - Le)m = with ui(xi, 0) = /(sin 2 ^/x~t) and 
(pt — L r )u r = with u r (x r , 0) = f(l — sin 2 



It is clear that the symmetry y — > 1 — y carries the left end to the right and vice- 
versa, so to simplify this discussion, we focus on the left end. We use x to denote 
xi, b to denote bo, and we choose a smooth cutoff function tp so that 



<p{x) 



'lforx e [0, (f) 2 -2t/] 
Oforx > (f) 2 -7/, 



where 77 > is small. With /i(x) = c(x)tp(x), we now focus attention on 

L = L bo + xh(x)d x . 
A solution to (fT5T ), pulled by via X£ satisfies: 



7T\ 2 



(a t -L)S = 0on[0,(^-J -277]. 

To build a parametrix for the heat kernel that has the correct boundary behavior 
near to y = 0, consider the initial value problem: 

{dt -L)u = with u(x, 0) = /(x) G C°([0, 00). (56) 

For our application / is obtained from / by pullback and multiplication by a 
smooth cut-off. Multiply the equation by k\(x, y) and integrate to obtain: 

t 00 00 

u M - 1 /^^(.H*,^ - /*.,,)/(»)*. 



or equivalently 

00 

= J k b t (x,y)f(y)dy. 
Q 

If 6 = 0, then as noted in Remark 13 we could replace k®_ s in the definition of A® 
with fe£l^, which shows that, as expected, u(0, t) = /(0) for all t > 0. 

It is straightforward to solve (1561 ) using the estimates from Section [8] Indeed, 
the solution can be expressed as a convergent Neumann series: 

00 

u(x,t) = {ld-A b t )- l g(x,t) = Y,(4) j 9, 

j=0 
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where 

g(x,t) = I k b t {x,y)f{y)dy. 



Let us denote the operator / — ► (Id — A^)^ 1 ^ by Q\f. Later on we shall need to 
distinguish the operator constructed near x = from the one near x = 1, and at 
that point we shall write them as Q t ° £ and Q h t l r , and denote their kernels by qfy , qfy, 
respectively. 

An immediate consequence of Propositions |9l and[His that if / G C*([0,oo)), 
then this sum converges uniformly in the topology of C°([0, T];C e ([0, oo)), for any 
T > 0, and Propositions [6l and [TOl show that for such data, 

lim i diu(-,t) = dlf for j < I, 

where the convergence is with respect to the C° topology. 

The regularity of these solutions shows that if / G C l ([0, oo)), then 

g € C°([0,oo);C £ ([0,oo)) nC°°([0,oo) x (0,oo)). (57) 

As a consequence of Proposition |7J we see that for < e < t, the sum defining 
Qtf actually converges uniformly in the C°° -topology. Hence 

u£ C ([0,oo);C^([0,oo)) nC°°([0,oo) x (0,oo)), (58) 

as well, and differentiating shows that 

(dt — Lb)u = xh(x)d x u <^=^ (dt — L)u = 0. 

We have now solved the problem 'exactly' near each endpoint, and the next 
step is to paste together these left and right solution operators to obtain an infinite 
order parametrix for dt — L, which has remainder term vanishing identically near 
the x = and x = 1. To accomplish this, define maps ^ and 4> r near x = and 
1 which put the leading part of L\vf mto tne model form xd%. Pulling back the 
operator in (fT8T > gives exactly the operator L in the intervals [0, |], and [J, 0]. Now 
choose cutoffs ipi,ip r ,ipo G C°°[0, 1] so that 

11 3 
ipi(x) = lforx G [0, — ], suppipi C [0, -], 
lb 4 

5 1 

(p r (x) = 1 for x G [— , 1], supper C [-, 1], 

3 5 
cp (x) = lforx G [0, -], supp^o C [0, -]. 
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and define the parametrix 



q t (x,y) = m{x)Qu (M x ), 4>i{y)) ( fi(y)\ ( l ) 'i(y)\+ 

(1 - ip {x))ql 1 r (4> r (x),(l) r (y))ip r {y)\(t>' r .{y)\. 

The heat kernel for L is determined symbolically to infinite order as t \ near the 
diagonal away from x, y = 0, or 1. Furthermore, the kernel and all its derivatives 
tend to zero like e~ c l l away from the diagonal and away from y = or 1. Hence 
in the overlap region, supp ipo n supp(l — <po), the two terms agree to all orders. 
Now set 

et(x, y) = [d t - (x(l - x)d 2 x + b(x)d x )}q t {x, y); 
by construction, this satisfies 

3 5 

suppe t (x,y) C [-, -] X [0,1] x [0,oo). 

Since this error term is obtained by applying derivatives in the first (x) variable to 
%(x,y),Corollary|5]showsthaty 1 - fco (l-y) 1 ~' ,1 e t G C°°([0, 1] x [0,1] x (0,oo)) 
and vanishes, along with all its derivatives, like e _C// * for some c > 0, on [0, 1] x 
[0, 1] x [0, oo) disjoint from the diagonal at t = 0. In a neighborhood including the 
diagonal it vanishes, in the C°°-topology, faster than 0(t N ) for any N G N. 
Given a function / G C m ([0, 1]), set 

i 

u Q (x,t) = / q t (x,y)f(y)dy. 



Propositions |9]and [JO] show that u (x,t) G C°([0, oo); C m ([0, 1])), and 

[d t - L]u (x,t) = v(x,t), lim u (x,t) = f(x), 

t->0+ 

where v(x, t) is smooth, vanishes in {[0, |] U [|, 1]} x [0, oo) and tends rapidly to 
zero as t — > + . This is even true for div(x, t) for all j G N. 

We now complete the construction of the solution operator for the generalized 
Wright-Fisher operator. Define, for e > 0, 

t-e 1 

Qt9= / Qt-s(x,y)g(y,s)dyds 



Et9 = ~j J e ts( x :y)9(y,s)dyds, 
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and write Q® simply as Qt. If 

g £ C°([0, oo) x [0, oo)) n C 2 ([0, oo) x (0, oo)), (60) 

then the estimates above imply that Q\g — > Q t g in the C 2 topology for any t > 0, 
and hence 

(8t - L)Q t g = lim (d t - L)Q\g = (Id -E t )g. 

The inverse of (Id —E t ) is an operator of the same type, and we write it as 
(Id —.Hi); here H t is represented by a kernel h t (x, y). Note that 

(Id -E t ) (Id -H t ) = (Id -H t ) (Id -E t ) = Id, 

so 

Ht = —Et + EtHt = —Et + HtEt ==^ Ht = —Et — E 2 + EtHtEt- 

The last identity shows that h t same regularity properties as et- In particular, it 
vanishes, along with all derivatives, to infinite order as t — > 0, and behaves like 
ybo-i along y = 0, and (1 — along y = 1. 

Setting g = (Id —H t )v, then ui = Q 4 g satisfies 

(dt — L)ui = v and limwifx-t) = 0. 

The true solution u is the difference 

u = u - ui = Qt(Id-E t + H t E t )(f (g> 5(t)), 

or equivalently, 



l t l 

u = 





1 I 1 1 

y Qt(x,y)f(y)dy + / y g 4 _ s (x,2;) ^ h s (z,y)f(y)dydzds. (61) 



This was derived under the assumption that g satisfies (l60l . If / G C°([0, 1]), 
then the fact that the solution can be represented the same way follows from the 
mapping properties of Q t . 

Definition 1. We let Q t f denote the operator on the right hand side in (|6TT >. and 
qt(x, y) its kernel. 
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Note that since y 1 ^ b °(l — y) 1_fcl /it decays rapidly decreasing as t \ in 
C°°([0, 1] x [0, 1]), the kernel q t already defines a complete parametrix for {dt — 
L)" 1 away from y = and y = 1. The contributions of h t along y = and y = 1 
are essential for Q t to satisfy the adjoint boundary conditions (fT3l , and hence to 
provide the solution operator for the adjoint problem. 

Putting together all the information we have obtained about the kernel for Q t , 
and using the uniqueness for the regular solution of (PT5T ). we can now state the 
fundamental 

Theorem 1. For each m £ N U {0} the operators Qt define positivity preserving 
semi-groups onC m ([0, 1]). The function u(x,t) = Qtf(x) satisfies (I15l >. More- 
over, for f £ C m , 

lim || Qtf - /||c™([o,i]) = 0. 

One of the reasons for having worked hard to establish convergence of the 
Neumann series used in the construction of the solution operator Q t is that we can 
estimate how good of an approximation the partial sums are. 

Proposition 12. Fix N > and define the operator Qt,N by pasting together the 
finite sums 

N 
j=0 

using the solution operators for the models at the left and right endpoints of [0, 1], 
and denote its kernel by qt^N- For any f G C m ([0, 1]), let u(x, t) denote the exact 
solution to (fT5T): for N > m and 2j < m, the function 

u N (x,t)= qt,N(x,y)f(y)dy (62) 
Jo 

satisfies 

||^[« J v(- J t)-«(-,t)]|| C o ([ o,i] ) <Ct j?± ^ 2i 

where the constant C can be estimated in terms of j, N, and the coefficients of L. 
In particular, when m = N = then 

sup |uo(^)*) — f( x )\ ^ C\ft. 
xe[o,i] 

We recall the changes of variables 

y/xi = sin" 1 s/x <fyi = sin" 1 Jy 



yfx~ r = sin 1 Vl — x ^fy~ r = sin 1 \J\ — y\ 
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differentiating we see that 



dyi _ sin yfij dy r 



sin 1 y/l — y 



(64) 



dy y/y(l-y) dy 



Vy(i-y) 



Substituting, we see that, for 60 an d b\ non-zero, the leading term in the solution 
operator qt is given by: 



Qt,o(x,y) = (fo{x)y e 1 ( — J e * ^ bo (x e y e /t 2 )ip e (y)—==L=+ 



10 The infinitesimal generator and long time asymptotics 
of solutions 

By TheoremQ] Qt defines a semi-group on C°; indeed, it also defines a semi-group 
on C m for every nonnegative integer m. To understand this C° semi-group better, 
and in particular to estimate the long-time asymptotics of solutions, we now seek 
a characterization of its infinitesimal generator A as an unbounded operator on C , 
including some features of its spectrum and a description of the behavior at x = 
and x = 1 of the elements in its domain. These facts will be proved using the 
various regularity results we have obtained. At various points in this discussion 
it will be necessary bring in the the adjoint operator I}, and in particular the in- 
finitesimal generator A* for the adjoint semi-group. Note that A* is an unbounded 
operator on [C°([0, 1])]', which we identify with A4([0, 1]), the space of finite Borel 
measures on [0, 1]. Section (T2] contains a more complete discussion of the adjoint 
semi-group. 

The first step is to note that the infinitesimal generator has a compact resolvent. 

Proposition 13. Let A be the infinitesimal generator associated to the C° semi- 
group defined by Qt- Then as an unbounded operator on C ([0, 1]), the spectrum 
of A lies in the left half-plane {A : Re A < 0}; furthermore, A has a compact 
resolvent. 

Proof. The first statement follows from the maximum principle, which implies 
that \\Qtf\\ c o < ||/||co. Next, observe that for any t > and / G C°([0,1]), 
Qtf £ C°°([0, 1]) The closed graph theorem and Arzela-Ascoli theorem now apply 
to show that Q t is a compact operator on C° for any t > 0, so the second statement 
follows from the results in Section 8.2 of Q. □ 





(l-ipo(x))y r I— J e ~tl; bl {x r y r /t )ip r {y) 
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Of course, the full characterization of A involves a detailed description of its 
domain. This will be based on a basic result from semi-group theory, due to Nelson: 



Proposition 14 (Nelson, Mill ). Let B be a Banach space and Z a closed operator 
on B generating a semi-group Tp If D C Dom(Z) is a subspace which is dense 
in B, and ifTfD C D,for every t > 0, then D is a core for Z, i.e. 



Z 



D- 



To apply this, recall that Q t C 2 ([0, 1]) C C 2 ([0, 1]) for t > 0, and also, if 
/ G C 2 ([0, 1]), then u = Q t f satisfies d t u = Lu for t > 0. 

Proposition 15. If A is the generator of the C° semi-group defined by Qt, then 



A - L le 2 ([o,l])- 

Proof. If / € C 2 ([0, 1]), then u(x,t) = Qtf(x) has one time derivative and two 
spatial derivatives, all of which are continuous on [0, 1] x [0, oo). Now integrate 
the equation satisfied by u to compute that 

Qtf-f i 



l r 

- / Lu(x, s)ds. 



t 

Since Lu G C°([0, 1] x [0, oo)) and equals Lf(x) at t = 0, we have 

^l—l = Lf + (l). 

This implies that C 2 ([0, 1]) C Dom(A) and on this subspace, Af = Lf. Since 
C 2 ([0, 1]) is dense is C°([0, 1]), the proposition now follows directly from Nelson's 
theorem. □ 

On the other hand, if / G Dom(A) then Lu G C°([0, 1]), so the one-dimensional 
version of "elliptic regularity" shows that / G C 2 ((0, 1)). In other words, the final 
characterization of Dom(A) involves only the description of its elements at the 
boundaries. 



The case where neither 6(0) nor 6(1) vanish 

It turns out that the results in case either 6(0) or 6(1) vanish are slightly more 
complicated to state, so for the moment let us suppose that < 6(0), —6(1). As 
before, denote 6(0) = 6q and —6(1) = 6i. 
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We begin by noting that there is a solution v to the adjoint equation L l v = 0, 
where 1} is given in (fTTT) . satisfying the adjoint boundary conditions ( fT3l >. Thus 

vq(x) = x^-^l ~ x) bx - l e B{ - x \ (65) 

where B[x) G C°°([0, 1]), and 

d x [x(l - x)vq](x) - b(x)v (x) = 0. 

The existence of i> follows using standard ODE techniques. 

Choose (p G C°°([0, 1]) with support in [0, 1) such that ip(x) = 1, for x in 
[0,i]- If /GC 2 , then 

l l 
(Lf(x))(p(x)v (x)dx = / f(x)L t (ipv )(x)dx. (66) 



Since C 2 is dense in Dom(^4), this identity also holds for the graph closure. If 
/ G Dom(A), and 5 is any small positive number, then 

i 

J [Lf(x)ip(x)v (x) - f{x)L\yvv){x)\ dx = -5(1 - 6)v (S)d x f(d). 

6 

Using (l66l ) and the asymptotic form of vq, we deduce that 

lim x b °d x f(x) = 0. (67) 

A similar argument using a cutoff function with support near to 1 shows that 

lim (1 -x) bl d x f(x) = 0. (68) 

x— +1 

Note that these are precisely the boundary conditions described in Section [3] They 
are also the ones described by Feller as defining a positivity preserving contraction 
semi-group on C°([0, 1]). 

If / G C°((0, 1)), then the definition of q t and the maximum principle imply 

l l 

lim / q t (x, y)f{y) dy = f(x) and lim / q t (x,y)f(x)dx = f(y). 

o o 

Using this and the sharp maximum principle for parabolic operators in one di- 
mension, see e.g. Theorem 2 in Chapter 3 of lfl2l . by a straightforward limiting 
argument we obtain a strict pointwise lower bound for q t '■ 
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Proposition 16. 7/" 6(0), and 6(1) are non-vanishing, then for t > ant/ x,y £ 
[0,1], 

q t (x,y)>0. (69) 

In this case, for each i > 0, the kernel q t (x,y) defines a strictly positive op- 
erator on C ([0, 1]). In other words, if / € C ([0, 1]) is non-negative and not 
identically zero, then Qtf(x) > for x G [0, 1] and t > 0. Consequently, we can 
now apply Theorem 23.1 from [ 10] (The Perron-Frobenius Theorem) to conclude 
that: 

Theorem 2. 7/6(0) ant/ 6(1) are non-zero, then the infinitesimal generator A is a 
compact operator on C°([0, 1]) with spectrum lying in {A : Re A < 0}. The only 
element in o~(A) on the imaginary axis is the point A = 0, and the only associated 
eigenfunctions are the constant functions. The function Vq, defined in (1651 ) spans 
the 0-eigenspace of A* . 

One consequence of the theorem above is that 

Ai = sup{Re(A) : A G a (A) \ {0}} < 0. 

Thus defining 

l l 
u(x,t) = J c\ t (x,y)f(y)dy, and c = J v (y)f(y)dy, 
o o 

where vq is normalized to have integral 1, then Theorem 2.1 in Section B-IV of ID 
implies the 

Corollary 6. Under the hypotheses of Theorem^ for each 5 G (0, |Ai|), there is a 
constant Ads > such that 

\\u(x,t)-c \\ c o <M s \\f\\ c oe- 5t . 
The case where either b = or b\ = 

We now turn to the characterization of Dom(A), and the corresponding decay re- 
sults for solutions, when either bo or b\ (or both) vanish. 

If b = 0, then by Proposition [BJ Af(0) = 0, for every / G Dom(i). This 
in turn implies that S(y) G Dom(A*), and that A*5(y) = 0. This is a new feature, 
since if 6o ^ 0, then clearly 5{y) G" Dom^*), and it is what complicates the 
discussion. Similar remarks apply at x = 1 if b\ = 0, of course. 
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More generally, if dfx G Dom(^4*) and A*d/j, = 0, then by elliptic regularity 
in the open interval (0, 1), the measure dfx has the representation g{y)dy where g 
satisfies L l g = and the adjoint boundary conditions, (fT3T >. 

d x [x(l - x)g(x)} - b(x)g(x) = 0, (70) 

so in particular g £ C°°((0, 1)). Writing b(x) = &o(l — x) — b\x + x(l — x)b{x), 
then solutions to (1701 have the form: 



e -B(x) 

9( x ) = C x l-fe (l_ x )l-fei 



(71) 



where C is a constant and 1? is a primitive of 6. Clearly, g(x)dx is a measure of 
finite total variation if and only if both bo, bi > 0. 

If both b = h = 0, then b(x) = x(l - x)b(x), where b G C°°([0,1]). The 
functions, 



Uo(X) 



x 

C J exp 
o 



b(z)dz 



L 



dy 



(72) 



are strictly monotonically increasing in (0, 1), and solve Lu = 0. Choosing C > 
appropriately, we can assume that 



u (0) = Oandiio(l) = 1. 
If bo = 0, but b\ > 0, then there is a non-negative solution no, such that 

uq(0) = 0, and lim (1 - x) bl d x u (x) ^ 0, 

while f 6o > 0, 6i = 0, then there is a non-negative solution uq(x) with 
uq(1) = 0, and 



(73) 



lim x °d x uo(x) ^ 0. 

E— >0+ 



If either bo or b\ vanish, then the solution operator to (fT5T ) can be expressed in 
a form analogous to (l34l . Suppose first that b vanishes at only one endpoint, say 
6(0) = 0, but 6(1) < 0. We can use the kernel ' to build a solution operator, q^, 
for the Dirichlet problem at x = with the regular boundary condition at x = 1: 

q t (x,y) = qP(x,y) + 5(y)co(x,t), 

where (dt — L)cq(x, t) = 0. As before, qf(x, y) > for x > 0, and moreover, 

i 



Co(x,t) = 1 - / q t {x,y)dy, 
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so Co > 0. A similar argument works if 6(0) > 0, but 6(1) = 0. Finally, if 
6(0) = 6(1) = 0, then we can write 

q t (x,y) = qf (x,y) + 5(y)c (x,t) + 5(1 -y)a(x,t), (74) 

where off (x, y) is positive on (0, 1) x (0, 1) and 

i 

c (x,t) + ci(x,t) = 1 - J q?(x,y)dy, 

- 

o 

In each of these cases, C°([0, 1]) splits into a finite dimensional subspace, in- 
variant under Q t , and an infinite dimensional complement, also invariant under Q t . 
For example, if 6(0) = 6(1) = 0, then 

C°([0, 1]) = C °([0, 1]) span{l, u }. (75) 

In all cases there is a corresponding splitting of the semi-group into the Dirichlet 
semi-group, Qff , and a semi-group on a finite dimensional space. The infinitesimal 
generator A D of Qf is the graph closure of L on the set of functions in C 2 ([0, 1]) 
which vanish at the appropriate end-point, or -points. The fact that qf is positive 
for < x, or x < 1, or < x < 1, respectively, implies as before that the 
semi-groups Qf are positive and irreducible. 

The next proposition follows from known results about positive, irreducible, 
compact semi-groups acting on Cq(X), where X = (0, 1), (0, 1], or [0, 1), see HI. 

Proposition 17. If either or both of the numbers 6(0), 6(1) vanish, then A D is 
compact. There is an element Ai G a(A D ) with \\ € (— oo, 0) and a unique cor- 
responding eigenfunction u\ which is smooth and positive in (0, 1). The remainder 
of the spectrum lies in Re A < Ai — rj, for some r] > 0. 

Proof. The compactness follows from the fact that Qf is compact for every t > 0. 
The existence of the eigenfunction u\ and the negativity of A i is obtained by ap- 
plying oscillation and comparison theorems for Sturm-Liouville operators. When 
6(0) = 6(1) = 0, comparison with operators of the form Md% + x ™_^ x ^ yields the 
existence of Ai < and a unique associated eigenfunction u\ £ Dom(A) with 
u\ > in (0, 1). The cases where only one of 6(0) or 6(1) vanish are somewhat 
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easier. If 6(0) = 0, 6(1) < 1, then the eigenvalue problem can written as 



d x [(l - xpe B{x) d x u} + -F -j-r- = 0, 

u(0) = and lim (1 - x) h d x u(x) = 0. 

Since (1 — x) bl ~ l is integrable near x = 1, we can apply standard oscillation 
theorems to obtain the desired conclusion. The result then follows from Corollary 
2.2 in Section B-IV of CD and the fact that A D has a compact resolvent. □ 

If 6(0) = 6(1) = 0, then any / G C° can be decomposed, according to (T73T ). as 

/ = /o + [/(0) + (/(l)-/(0))«oL 

and then 

Qtf = Gf/o + If (0) + (/(l) - /(0))«o]- (76) 

The proposition above and Theorem 2.1, in Section B-IV of [1J imply that there is 
a constant, independent of / so that 

||Q f D /o|| c o<Ce A ^||/|| c o. 

Comparing the representations of Q t f in (1741 and (1761 ) shows that 

c (x, t) = l- u (x) + 0(e Alt ) and Cl (x, t) = u (x) + 0(e Xlt ). (77) 

There are similar results for the other two cases, where only one of 6(0) or 6(1) 
vanishes. In this case the solution tends asymptotically to the constant /(0), if 
6(0) = 0, and /(l), if 6(1) = 0. 

11 The Resolvent of A 

The Hille-Yosida theorem states that if A is the infinitesimal generator of a con- 
traction semi-group Qt, then the right half plane belongs to the resolvent set of 
A. For A with positive real part the resolvent, (A — A)^ 1 , is given by the Laplace 
transform of Q t : 

oo 

(A - AY 1 = J e~ xt Q t dt. (78) 
o 

If A is the C°-graph closure of a generalized Wright-Fisher operator, L, then Q t , 
defined above, is a contraction semi-group on C°([0, 1]), and therefore the right 
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half-plane is in p(A). In this section we consider the higher order regularity of 
solutions to (A — A)w = f. 

Our regularity results show that if / G C m ([0, 1]) then the solution, u to ([T5l 
satisfies 

uGC°([0,oo);C m ([0,l]))nC oo ([0,l] x (0,oo)). (79) 

We can therefore differentiate the equation satisfied by u with respect to x to obtain 
that, for 1 < j < m, and t > 0, 

d t [diu] = x(i - x)d 2 x [diu] + (b( x ) + j(i - 2x))d x [diu] + Cj {x)[diu\, (80) 

for a function Cj G C°°([0, 1]). The operator 

L {j] = x(l - x)d 2 x + (b{x) + i(l - 2x))8 x , (81) 

is a generalized Wright-Fisher operator. Applying a standard extension of the max- 
imum principle to equations with a zero order term, see Theorem 4 in Chapter 3.3 
of lfT2l . we easily obtain 

Proposition 18. For m G N, there are constants C m ,^ m > 0, so that for f G 
C m ([0, 1]), the solution to (1151 ) satisfies 

IM,*)llc»([o,i]) < C m e^*||/|| cm([0il]) . (82) 

Combined with the continuity result in Theorem [T] this shows that Q t defines a 
semi-group on C m ([0, 1]) with || Q t ||c m < C m e Mm *. We let the infinitesimal gener- 
ators be denoted by A m . By Nelson's theorem, these can be taken as the C m -graph- 
norm closure of L acting on C°° ([0, 1]). As Q t is a compact operator on C m ([0, 1]) 
for t > 0, the operators A m are compact for all m G N. 

Evidently the resolvent set of A m contains the half -plane {Re A > /i m }. Be- 
cause A m is a compact operator, if A G a(A m ), then there is an eigenvector 
u G Dom(^4 m ) so that 

Lu = An. (83) 

Because Q t u = e xt u G C°°([0,1]), we see that the eigenvectors all belong to 
C°°([0, 1]), and therefore an eigenvector of A m is also an eigenvector of A, and 
vice versa. This shows that, as a point-set, 

a(A m ) = a(A) C {A : Re A < 0}. (84) 

Theorem 3. If m G N, then for f G C m ([0, 1]) and A, with positive real part, the 
solution w G Dom(A) to the equation 

(A - L)w = f, (85) 
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belongs to C m ([0, 1]) and ifReX > p m , then 



\\w\\c m < 7T r. (86) 

I A ~~ Mm I 

Proof. The Hille-Yosida theorem shows that the estimate (l82l) implies that {A : 
Re A > Hm} belongs to the resolvent set of A m . For / G C° we set 



-Af 



Qtfdt. (87) 



This is an analytic C°-valued function in {A : Re A > 0}. If / G C m , then in 
{A : Re A > p m }, this equals (A — A m )~ l f, which is an analytic C m -valued 
function in, p(A m ), the resolvent set of A m . As noted in (l84l . p(A m ), includes 
{A : Re A > 0}, which completes the proof of the theorem. □ 

This is not, in any real sense, an elliptic estimate, as it only shows that the 
solution w is at least as regular as /. Of course in the interior of the interval, w has 
two more derivatives than /, but this may not be true, in a uniform sense, up to the 
boundary. 



12 The adjoint semi-group 

It is a consequence of the boundary behavior of qt, which follows from Corollary \5\ 
that for t > 0, and each x £ [0, 1], 

L*q t (av) G ^([0,1]) (88) 
The results in Section |9] show that, if / G C m ([0, 1]) for m > 0, then 



c(x,t) = J q t {x,y)f(y)dyeC°([0,oo) t ;C m ([0,l} x )) 



i 



and satisfies dtu = Lu with u(x, 0) = f(x). Thus for / G C 2 ([0, 1]), uniqueness 
for this initial value problem shows that, for t > 0, 



Lu(x,t) = J q t (x,y)Lf(y)dy. 



46 



By Corollary [51 qt satisfies the adjoint boundary conditions (fT3T >. so the integrabil- 
ity of Ly q t , implies that we can integrate by parts to conclude that 



Lu(x,t) = J L y q t (x,y)f(y) 



dy. 



Since 



d t u = J d t q t (x,y)f(y)dy, 
o 

it follows immediately that, for all / G C 2 ([0, 1]) and t > 0, 

t 

f (d t -L t y )q t (x,y)f(y)dy = 0, 
o 

and hence we conclude by a straightforward limiting argument the following: 
Theorem 4. If g e C° ([0,1]), then, fort > 

i 

v{y,t) = j q t (x,y)g(x)dx, 
o 

satisfies the boundary conditions (113I ). and solves the initial value problem 
(d t - L l )v{y, t) = and lim v(y, t) = g(y). 



(89) 



The dual space of C°([0, 1]) is naturally identified with M([0, 1]), the space 
of Borel measures with finite total variation on [0,1]. The dual semi-group, Q' t , is 
thus canonically defined on this space by 

(Qtf,df J ) = (f, Q! t dii). 

However, C°([0, 1]) is not a reflexive Banach space, so the dual semi-group is 
weak* -continuous, but not necessarily strongly continuous. The infinitesimal gen- 
erator A of Q t has a canonically defined adjoint, A*, whose domain is defined by 
the prescription: a measure dv E Dom(>l*) if there exists a constant C so that for 
every / £ Y)om(A), 



Af{x)dv{x) 



<C||/|| c o. 
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The subtlety is that Dom(A*) may not be dense in ,M([0, 1]). Following Phillips, 
see Q, we define the adjoint semi-group as 

Qf = Q' t \ M<3j where M & = Dom(A*)nA* Dom(A*). 

Phillips shows that Qf is a strongly continuous semi-group on A4 & , with infinites- 
imal generator: 

A Q = A* lDom(A*)ncM© • 

Thus our task is to identify M. G . 

By TheoremHl measures of the form d[i = g{y)dy with g £ C 2 ((0, 1)) are in 
Dom(yl ) provided L l g G L 1 and g satisfies the adjoint boundary conditions ( PT3T ). 
The closure of all such measures with respect to the topology of A4([0, 1]) is 
L 1 ([0, 1]). Any classical eigenvector of L l satisfying (fT3T > belongs to El- 
liptic regularity implies that a distributional solution to L l d\i = dv is absolutely 
continuous on (0, 1), and hence 

^([0, 1]) CM G . 

In fact, except for the possibility of atomic measures at and/or 1, these two spaces 
are equal. 

If 6(0) ^ 0, then for any / E C 2 ([0, 1]), 

(Lf,S(y)) =6(0)^/(0). 

Since the right hand side does not represent a bounded functional on C°([0, 1]), we 
obtain that 5 ^ Dom(^4*). A similar calculation can be done at x = 1. On the other 
hand, if 6 vanishes at either end, then: 

1. If 6(0) = 6(1) = 0, then the nullspace of A & is spanned by 5(y) and 6(1— y), 
both of which belong to .M . 

2. If one of 6(0) or 6(1) is non-zero, then ^4® has a 1-dimensional nullspace 
spanned by 5(y), if 6(0) = 0, or 5(1 — y), if 6(1) = 0, and this nullspace 
again lies in M® . 

Summarizing these observations we have proved the 

Proposition 19. • 7/6(0) and 6(1) are non-vanishing, then M® = -^ 1 ([0, 1]). 

• //6(0) = 6(1) = 0, then M G = ^([0, 1]) span{%), 5(1 - y)}; 

• 7/6(0) = 0, 6(1) ^ 0, then M G = L^fO, 1]) span{<5(y)}; 

• i/6(0) / 0, 6(1) = 0, then M G = ^([0, 1]) span{5(l - y)}. 
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The spectrum of the operator ^4° equals that of A, see Theorem 14.3.3 in Q; 
therefore we have the same sort of asymptotics for solutions to 



Theorem 5. If < 6(0) and 6(1) < 0, then there exists a X\ < 0, so that the 
solution v(x, t) to d89l . with initial data f E -^([O, 1]) satisfies 



v(x,t) = c v (x) + 0(e Xlt ), where c = J f(y) 



dy. 



Here vq E Dom(A*) is the positive solution to L^vo = 0, normalized to have 
integral 1. If 6(0) = 6(1) = 0, then there is a X\ < such that 

v(x,t) = c Q d(y) + d$(l - y) + 0(e Al *), 

l l 

w/iere c = J (1 - u (y))f(y), and c\ = j u {y)f(y) dy. 
o o 

where uq is defined in (1721 ) awif d73l . Finally, if only one ofb(0) or 6(1) vanishes, 
then 

fc5(y) + 0(e^), 6(0) = 0, / n , n , 

^' t) = \c5 ( l- y) + ( e^), 6(1) = 0, C = J U ^^ d y- 

Proof. The first case is immediate from the fact that is an isolated point in the 
spectrum of ^4®, and all other eigenvalues have strictly negative real part. The 
second case follows from dTTT ) and the analogous fact about the spectrum of A D * ; 
the last two assertions are similar. □ 

As a simple special case of our regularity theorem for solutions of the back- 
wards equation we have: 

Corollary 7. lfu{x, t) is a solution to 

d t u = x(l - x)d 2 x u u(x, 0) = f(x), (90) 

for f E C 2l ([0, 1]), then, for (j + k) < I, the functions &[L k u are continuous on 
[0,1] x [0,oo). * 

If v solves the forward Kolmogorov equation: 

d t v = d 2 y [y{l - y)v] wifh«(y,0) = g{y) eC 2i ([0,l]), (91) 
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then it is a simple calculation to see that u(x,t) = x(l — x)v(x,t), solves the 
backwards equations, with u(x, 0) = x(l — x)g(x), and u(0, t) = u(l, t) = 0. 

The theorem shows that x(l — x)v(x, t) is therefore a C 21 function on [0, 1] x 
[0, oo), vanishing at and 1. This easily implies that v itself is in C 2l ([0, 1] x 
[0, oo)). If we let Ly^ F g = d 2 x(l — x)g, then this regularity shows that for any 
1 < j < I 

(92) 



a r'ti" 1 ),, _ r*i „ 



We can integrate the forward equation, and use this formula, repeatedly integrating 
by parts, to obtain 



v(x,t)=g(x) + J L WF v(x,s)ds 



= g{x) + {s- tfLyjpvix, s) 



s=t 



r 

+ J (t- s)L$ F v(x,s)ds (93) 



j=0 



(1-1)1 



s=0 o 

(t - s) l L\k F v(x,s)ds. 



Using the same argument, we can show that if / G C 2Z ([0, 1]) and u solves 
the initial value problem: (pt — L)u = 0, u(x,0) = f(x), for L a generalized 
Wright-Fisher operator, then 



l-l 



u(x 



N ^ PVf 1 



(t - s) l L l u(x, s)ds. 



3=0 o 

Applying Proposition [18] show that there is a constant M 2 ; so that 

\L l u(x,s)\ < M 2 ie ms 

and therefore: 



«(a;,t) - Yj 



j=0 



M 2l \\f\y ie ^H l 
i\ 



< 



(94) 



(95) 



(96) 



A Appendix 

Proof of Proposition [9j 
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Consider the functions 

oo 


defined when t > s; these satisfy 

(d t — Lb)v = and v (x, s; s) = xh(x)d x g(x, s). 
By the maximum principle and the derivative bound on g it follows that 

\v(x,t;s)\ < suv\xh{x)M / yfsx\ < \/LM\\h\\oo, 

and hence 

t 

\A h t g(x)\ < J M^L\\h\\ 00 ^= = 2M^Tt\\h\\ 00 . 
o 

This establishes (1501) for j = 1, with do = 1. Now apply Lemma[5]to get 

C 6 ||u(-,s;s)||oo 



\d x v(x,t;s)\ < 



y/x(t-s) 



so that 



\d x A b t g(x)\ < 



1 CtMVZWhWocds 



\Jxs(t — s) 

u 

CwMy/LWhWn 



setting d\ = tt, this is (IBTI ) for j = 1. 

Assume that we have chosen do, . . . , dj—i. Now write 



/•oo 

Vj(x,t;s)= / k b t _ s {x,y)ydyv j -i{y,s)dy. 
Jo 



Using d5TT > for j — 1, this is bounded by 

f'°° -i 
d^MCl-WIWhW^y / kl s (x,y)s^dy 



Id^MCi-^y/LllhlUHt. 
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since J °° k b _ s (x, y) dy = 1 for s < t, which establishes (1501 ) for j. 

Using the induction hypothesis one more time, insert the result into the estimate 
of Lemma[5] Noting that l/(t — s + y/x(t — s)) < 1/ y/x[t — s), we get 



J Jxit — s) 

v v / 

_ dj^MiCbVLWhWooyt 1 ^ [a^rda 



ft i vr 



Evaluating the integral shows that we should set 



dj = \pHdj—\ 



and this proves (151b for j, thereby completing the induction. A straightforward cal- 
culation shows that dj is given by the formula in the statement of the Proposition. 

□ 

Proof of Proposition [70| 

We start with j = 1. These results follow from d48l ), the maximum principle, 
and Proposition |9l Combining these ingredients shows that when t < T, 

\d e x A b t g(x)\<V~t[2 e+1 \\xh\\ ce ](M + \\g\\ 

and, 

\di +1 A b t g{x)\ < 



d 1 MC b+t VL\\,^ + ^ + i ||lft||cw(1 + ^(ll^ll^ + M), 



x 



which establishes (1531 ) and (1541 ) when j=l. The main issue is to see how Dj de- 
creases as j increases. We assume that these estimates have been established for 

Applying (1481 ), we see that 



dUA b t y 9 (x)=AW(di(A b y- 1 g)(x)+ 

g t oo 



P =i 







52 



Using the induction hypothesis and the maximum principle, the first term on the 
right here is bounded by 

2£>j_ivT||/t|| 00 t2 , ._j 

^ ( M + \\9\\c^)[ C T,L,£,b\\ h \y\ ■ 

The terms in the sum with 2 < p < £ can be bounded using the maximum prin- 
ciple and Lemma [6l while the term with p = 1 is controlled using the induction 
hypothesis. Thus altogether, the sum is bounded by 



(2 £ \\xh\yy\\g\yp 2D j _ 2 P+ 1 2{M + ||9llc^)ll^lloo[C^ L/i J/i|| c 



J'! i(j - 2) 

So long as 



1 2Dj_ 2 



D ^>^\j V -±f)^ (97) 



then there is some C' T L t b so that (I53T ) holds for all j. That (I97T ) holds follows 
easily from the proof of (l54l ). to which we now turn. 
We next apply (|48T ) to see that 



t oo 



di +l {A b t y 9 {x)=d x J J k b + e s (x,z)zh(z)di +1 (A b s y- 1 g)(z)dzds+ 



t oo 

tyd x J J k b + i s (x,z)d y (zh(z))di(A b s y- 1 g)(z)dzds+ 



f t oo 



dzds. 



o o 



By the induction hypothesis and LemmaO the first term is estimated by 



J C b+ i\^L\\h\\ 



j 2 



rn + ^ T & ( M + \\9\y,~)(CT,L,t,b\\h\\ C ' + i) j - 1 
= \/TrJJj-it 2 — — 

3 r(i±±) ^ 

and the second term by: 

yPHiD^th r(|) (m + WgWc^ycT^WhWc^y- 1 
j-2 r(i±i) Vi 
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Using the induction hypothesis, Lemma[6]and the maximum principle, the last term 
is bounded by 

( M+ Mc^)(c^ LAb \\h\y +1 y . 

Hence if Da satisfies the recursion relationship in the statement of the proposition, 
then Dj-i satisfies d97l) , for j > 2 and there exist constants CT,L,e,b and C' T L ^ b 
so that 423) and §4$ hold for all j > 2. ' □ 
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